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Abstract. We consider a strongly interacting Bose-Einstein condensate in a spherical 
harmonic trap. The system is treated by applying a slave-boson representation for 
hard-core bosons. A renormalized Gross-Pitaevskii theory is derived for the condensate 
wave function that describes the dilute regime (like the conventional Gross-Pitaevskii 
theory) as well as the dense regime. We calculate the condensate density of a rotating 
condensate for both the vortex-free condensate and the condensate with a single vortex 
and determine the critical angular velocity for the formation of a stable vortex in a 
rotating trap. 
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1. Introduction 

In this paper we sliall study two aspects of a strongly interacting Bose gas at high 
density. One is related to a consistent treatment of a strongly interacting Bose gas in 
terms of an effective Gross-Pitaevskii (GP) Eq. with renormalized parameters. The 
second aspect is related to the formation of a vortex in a trapped condensate in the 
presence of strong interaction. 

The stationary form of the conventional GP Eq. 



describes the condensate order parameter $ of a Bose gas in a trapping potential at zero 
temperature, where /i is the chemical potential and g the repulsive coupling constant 
[H [21 [3] . In the absence of a trapping potential, a solution of Eq. ([1]) is given by 



This describes a linearly increasing condensate density (the latter is proportional to 
with respect to the chemical potential. Although it takes the repulsion into account by 
a factor 1/g which is decreasing with increasing coupling constant g, the saturation 
of no cannot be seen in this solution. From the physical point of view, in a realistic 
description for large densities, the particle density must saturate because there is a finite 
scattering volume around each particle. Furthermore, for increasing particle density, 
the condensate density should reach a maximum and for even larger densities, decrease 
again until its total destruction, because of the increasing interparticle interaction. This 
behaviour has also been found by variational perturbation theory [1] and diffusion Monte 
Carlo calculations [5]. In other words, the strong effect of the repulsion in a dense 
condensate is not really described by the conventional GP Eq. In order to describe 
condensates at higher densities, the second order term in the low-density expansion of 
the energy density has been taken into account which lead to a modified GP theory 



Although in many experimentally realized situations the BEG is in the weakly 
interacting regime where it is well described by GP theory, it might be possible to 
reach the strongly interacting regime. The main problem at high particle densities 
is the instability of the Bose gas by the formation of molecules due to three-particle 
interactions [8]. Here, we will assume that molecule formation does not occur. This 
might be unrealistic for some systems, but in others it is not, e.g. for Bose gases in 
optical lattices. 

It has been shown that the slave-boson approach to a hard-core Bose gas provides 
a mean-field Eq., similar to the GP Eq., that leads to a saturation of the condensate 
'^o ^ 1 01 Uni [in Il2l Us] . Here we will discuss a simplified version of this mean-field 
equation which is the same type of nonlinear Schrodinger Eq. like in the GP approach. 
However, in contrast to the latter the parameters are renormalized such that Eq. ([2]) 




(1) 




(2) 



[IlEl[6l[7]. 
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becomes 

I* |2 - 

9r 

where the renormahzed parameters /xr and (/r are functions of the bare parameters /i 
and no- At the phase transition between the non-condensed Bose gas and the BEC /iR 
vanishes, then increases hnearly with increasing /i, reaches a maximum and decreases 
again until the condensate is destroyed totally due to strong interaction effects. We use 
this approach to calculate the condensate density profile of a trapped BEC. 

The formation of vortices in a rotating condensate is understood as a local 
destruction of the condensate. From the depletion effect due to strong interaction it 
can be anticipated that the tendency of vortex formation is enhanced in a dense Bose 
gas. This implies a reduction of the critical angular velocity Qc with increasing density 
by the interaction. Vortex formation in rotating traps in the region beyond the validity 
of GP has also been studied by the modified GP theory and by variational Monte Carlo 
methods [7]. 

The paper is organized as follows: In the following section we survey the results of 
the slave-boson approach for a hard-core Bose gas on a lattice and derive a corresponding 
mean field theory for a continuous system in an external trap potential. Then we derive 
the renormalized GP Eq. and evaluate solutions with and without a straight vortex and 
determine the critical angular velocity for the vortex formation. 



2. The Model 



The slave-boson representation of a strongly interacting Bose gas is based on the idea 
that the bosons fill the space with finite density, where each particle occupies a lattice 
cell related to a sphere with radius (s-wave scattering length) [10]. Singly occupied 
and empty sites are described by two complex fields and e^, respectively [H]. Here, 
b* and b^ represent a creation and annihilation process of a boson at site x, while e* 
and Cx represent the creation and annililation of a "hole". In the functional integral 
representation [H], the grand canonical partition function in classical approximation is 
given as 

exp -p bleJx,x'el,bx' + ^^/i^^lfexP 

x,x' X 

X W\ex? + \bx? - I)d6,d6:de,de: , (3) 

X 

where 1//? = k-^T is the thermal energy. In a (i-dimensional lattice, 

J —J/ {2d) if X, x' nearest neighbours 



tx,x' 



else 



is the nearest-neighbour tunneling rate and fix = fJ' — Vx is the space-dependent chemical 
potential that includes the trapping potential Vx- Instead of fix we will write only fi 
subsequently and assume implicitly that the effective chemical potential can depend on 
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space. The fields e and h are dimensionless. Moreover, we rescale all physical energies 
by a multiple of the hopping rate a J to obtain dimensionless quantities: 



i ^ i' 



-t 



(3^f3' = aJ(3 



(4) 



aJ aJ 

The hopping term is quartic in the field variables, due to the fact that a hopping process 
is characterized as an "exchange process" of a boson and an empty site. The 5 function 
enforces the constraint that each lattice site is either singly occupied or empty but 
excludes a multiple occupation. 

A similar approach has been applied to the Bose-Hubbard model which allows 
multiple occupation, by introducing additional fields, one for each occupation number 
p!2| [T3] . The well-known zero-temperature phase diagram with Mott-insulating phases 
for integer lattice fillings = 0, 1, 2, . . . ("lobes") was found. In contrary, our hard-core 
Boson model is restricted to the two lattice fillings n = 0, 1 which is a simplification, but 
contains all relevant aspects of Bose-Einstein condensation with repulsive interaction. 
A second simplification of our approach is, that we neglect quantum fluctuations by 
treating the grand-canonical partition function in classical approximation. This allows 
us to integrate out the constraint exactly. On the other hand, it restricts the applicability 
of our approach to non-zero temperatures. 

Two new fields are introduced by a Hubbard-Stratonovich transformation, the 
complex field $a; which describes the condensate wave function, and the real field (p^ 
which is related to the total density of bosons [TO] : 



Z = / exp I - [ E (l - i') Z, + 11^1 



2ip^ + 1 $ 




Y[5{\b.? + |e.r - I)d6,d6:de,de:d<l>,d$:d¥., . 



(5) 



Integration over the fields $ and leads back to Eq. ([3]). On the other hand, the fields 
hx and can be integrated out in ([5]) because they appear in the exponent as quadratic 
forms. This leads to the partition function 



with the kinetic part of the action 

x,x' 

and the potential part 
/ 



(7) 



1 / / 
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The condensate density can be identified with 



no 



(1 + l/aY ' 

hke argued in Appendix A, and the total particle density is given by the expectation 
value Ho] 



ntot 



{"fx) + ^ • 



(10) 



We apply a saddle-point approximation to the integration in Eq. ([6]). This is controlled 
by the minimized action, which means that we have to solve the equation 
dS 



. 



This yields the mean-field equation 



i' 



d 



,5(1$. 



log^a; 



. 



In order to derive a mean-field equation which is applicable to a continuous trapping 
potential, we perform the continuum approximation of flTT]) . If the field $r is varying 
only very slowly between neighbouring lattice sites, we can approximate 



and perform the substitution 

y^(^ 5x,x' + tx,x') 

x' 

This leads to the equation 
Jo" 



1 



l/a 



Sx.x' + 



1 + 1/a 



+ ix.x' 



6 



a)J — 



{1 + l/af 



d 



13 d{\^r) 



■logZ(r) 



$(r) = 



(12) 



with 



Ziv) 



sinh 



/3V(¥^r + /xV2)2 + |$(r)P 



dv9r 



_/3y(y,r + /iV2)^ + |$(r)P_ 

for the spacially dependent order parameter $(r) in a three-dimensional space and a 
is the lattice constant of the discrete system. In the continuum a looses its identity 
as lattice constant, but describes a characteristic length scale in Eq. (fT2ll . and can be 
interpreted as the spacial extension of a boson. Thus, it should be of the same order of 
magnitude as the s-wave scattering length [Ij. Eq. f|T2|) is the analogue of the GP 
Eq. in the case of our slave-boson approach. The parameters can be identified with 
those of the conventional GP Eq.: The mass m of the particles is given by the hopping 
constant J and the lattice constant a via 



2m 



6 



(13) 
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3. Renormalized Gross-Pitaevskii Equation 



The continuum limit of the action defined by Eq. ([71) and ([H]) is 



S 



/5'$*(r) 



a 



(1 + 6 



1 + 1/a 



$(r) -logZ(r) \ d'^r . (14) 



Applying the variational principle 
dS 

, 



we obtain the full mean-field Eq. (fT2l) directly. If the order parameter $ is small, we 
can expand the potential part of the action up to fourth order: 

<l>p-ilogZ(r) =ao(/i')+a2(/i')l<^'r+J«4(/i')l<^'r+(^(l<^'r) ,(15) 



1 + 1/a 



where we have introduced the coefficients 



- ^log^(r) 



$=0 



d 



(3'd\^ 
1 



2 



logZ(r) 



/3' (a|$|2)2 

Further, we introduce the rescaled field 

1 



<s>=o 
log Z{r 



1 + 1/a 



(16) 
(17) 
(18) 



<I'=0 



$R(r) 



<|.(r). 



1 + 1/a 

which we identify with the condensate wave function of our renormalized GP theory. 
We now introduce the renormalized coefficients 

(1 + aY 



and 



AiR(/^, J) 



gnifJ', J) 



a 



[1 + a) 



a- 



- Ja2(/i') 



-Jai{ij!) 



After neglecting the term of order |$|^ in the expansion ( fTSl) . we get the renormalized 
Gross-Pitaevskii (RGP) Eq. 

Jo" 



6 



-V - /iR(/i, J) + ^R(/i, J)|$R(r) 



$R(r) = . 



(19) 



It has the same form as the conventional GP Eq. ([T]), when /i and g are replaced by 
and gR. Their /^-dependence is plotted in Fig. [H In the case of a trapping potential, 
where the chemical potential /i is space-dependent, /xr and (7r are space-dependent 
aswell. While ^^r is always positive, /ir can change sign. A BEG exists if /xr > 0. The 
phase transition between the BEG and the non-condensate phase, i.e. the point at which 
the condensate order parameter vanishes, is given by the relation /iR(/i, J) = 0. 
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4. Results 



4.I. Zero-temperature result 



In the zero-temperature limit we can integrate out the (/9-field in Eq. (8) exactly by a 
saddle-point integration, as shall be shown in this paragraph. Therefore, we write 

1 



where 



and 



-/3'/±(¥',l*P) 

dip 



(<^ + ^)2 + |$|2 



It is possible to perform a saddle-point approximation by expanding the functions f± 
up to second order in (p around their minimum (po: 

In the limit of large /3, the saddle-point integration becomes exact and yields 



z. 



± 



2 9(^2 



The minimum is found to satisfy the equation 

2/ Wo 



and we have 

/±(^o) = ^0 



2 ^ 

We find the following zero-temperature result for a translational invariant condensate 
{^x = $ = const) from the mean-field Eq. (JT2l) : For the condensate density we find 

° (1 + \ else , ' 

and the total particle density given by Eq. ffTOl) is 

if /i < - J 

ntot = ^o + l = \ Ul-^) if - J < /i < J (21) 



2 V 

1 if J < /i . 

The solution of the mean-field equation is plotted in Fig. [2] for zero-temperature and 
near the critical temperature Tc where the BEC breaks down. In the dilute regime at 
T = 0, the chemical potential can be written as = — J + A/x, where A/x <C J. In this 
limiting case we find /zr = A/i -|- O^Afi"^) and g^ = 2J, which is consistent with the 
conventional GP Eq. with a shifted chemical potential. 
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4-2. Vortex-free trapped condensate 

Assuming a dense condensate, where the repulsive interaction between bosons dominates 
their kinetic energy, we neglect the differential term in Eqs. ( fT2i) and ( |T4l) . This is called 
the Thomas- Fermi (TF) approximation [1] . In the following we use a spherical trapping 
potential 

Vir) = "^uiy . (22) 



In typical experiments, the oscillator length dho = \J'h/mujYio is of the order of a few //m 
[1], where Who is the trap frequency measured in Hz. Considering, for instance, ^^Rb 
atoms near a Feshbach resonance [15], we can study a Bose gas in a dense regime with 
a scattering length ~ a ~ 200nm. In our calculations we choose the parameters 

/5' = 1, = 36.93, ^ = 0.1215, (23) 

AlCUho "ho 

and keep the hopping constant J fixed. Thus all energies can be scaled with J. 

To calculate the profile of the condensate density in a BEC without vortex, we solve 
the TF equation 

(1 + l/a) - (1 + l/«)2^^^ log Z{v) = . (24) 
In the RGP approximation f|T9l) . the solution is 

=5I="(l + l/a)W)- 
Solutions for typical values of the parameters are plotted in Fig. [31 The results we 
get from the renormalised GP approximation shows only small deviations from the 
numerical solutions of Eq. (IMl) . We find a condensate depletion at the trap center 
for /i' = 1. This is due to the fact that the condensate is partly suppressed by strong 
interaction effects [U [5l [H]. For fj,' = 2 the condensate is completely destroyed at the 
trap center. We find particle numbers in the condensate of the order A^^o ~ 10'^ . . . 10^. 

We note that the total particle density ntot is much larger than the condensate 
density no and takes values of about 0.5 at the trap center. Thus the interaction 
between the non-condensed and the condensed part of the Bose gas plays a significant 
role. This implies that the conventional GP Equation, which neglects the non-condensed 
part, is not reliable in this parameter regime. 



4-3. Rotating condensate with a single vortex 

In the case of a trap rotating about the 2;-axis with an angular velocity Q, one must 
include the additional angular momentum term —QLz^{r) to the left hand side of 
the differential equation ( |T2i) . where is the ^-component of the angular momentum 
operator. This term must also be kept in the TF approximation. The condensate 
wave function may develop a vortex then. We assume here a straight single vortex 
along the 2;-axis. This can be described by using cylindrical coordinates and the ansatz 
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$(r) = 0(r_|_, z)e*'^, where r± is the distance from the 2;-axis and (f the polar angle. The 
angular momentum operator is given as = This gives rise to an additional 

term {Vt/aJ - / a{Qr\))\^{v)\'^ / {I + l/af in the action Instead of Eq. (El) we 
have to solve 

+ ^/") + - If) - <i + ^'"-^'imyw '"^ = ° ■ 

and the solution in the RGP approximation is 

1 



^ ) + (1 + l/afa2{^i') 



a6r , aJ , 



■ (27) 



(1 + l/a)%4(/i'; 

A condensate that is rotating with given angular velocity Q forms a stable vortex, 
if its total energy is lower than that of a vortex free condensate. This is equivalent to 
the condition 

S^'^'^iQ) -S <0 (28) 

which can be checked numerically by using the TF approximation of Eq. f|T^ for a 
condensate with vortex 5'pT*(^^) and without vortex Stf- The critical angular velocity 
Qc above which the vortex is stable is plotted against the number of condensed bosons 
A^o in Fig. HJ where Nq is given as 

No = I ^I^RWpd^r. (29) 

The RGP approximation is in good agreement with the results of the full mean-field 
equation. The decreasing critical angular velocity for higher values of A^^o indicates that a 
high interaction energy favours the formation of a vortex. It agrees with results derived 
from the CP Eq. by perturbation theory [16] as well as numerically [T7] . 

Typical solutions for shapes of condensate density profiles of BECs with a stable 
single vortex are shown in Fig. [51 In contrast to the case without a vortex, the 
condensate is always completely destroyed at the trap center, a feature that also shows 
up in the conventional GP approximation [T] . Again, the RGP approximation is in good 
agreement with the numerical results from (l26l) . 



5. Conclusion 



The slave-boson approach allowed us to study the condensation of a trapped high- 
density Bose gas in a regime where the conventional Gross-Pitaevskii approach is not 
valid. Starting from the saddle-point approximation, we have derived a renormalized 
Gross-Pitaevskii Eq. with a space dependent coupling constant. This provides good 
results in comparison with the more complicated slave-boson saddle-point calculations. 
At high densities, we have found a depletion of the condensate at the trap center due 
to the interaction between the condensate and the non-condensate part of the Bose gas. 
This feature is not covered by the conventional Gross-Pitaevskii Eq. The behaviour 
of the critical angular velocity for the formation of a single vortex agrees qualitatively 
with previous results in the literature but supports also the formation of a vortex for 
increasing Nq. 
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Appendix A 



For a Bose system with creation operators and annihilation operators at a lattice 
site appropriate definition of the condensate density is [3] 



no :-- 



lim 



X— x'— ►oo 



ax' 



(30) 



In our slave-boson representation, a creation process of a particle is associated to the 
product 6*62, and an annihilation process to 6*6^., thus 



no 



lim (6*.e^e*,6^/) 



x—x' ^00 



(31) 



Here, the expectation value is defined with respect the functional integral ([5]) by 

(■■■) = 1/ ■■■exp[...]n5(l&.f+ |e.f -l)dM6:de,de:d$,d$:d</.,. . 

We are interested in the connection between the correlation function and the 

condensate density. For this purpose we integrate out the field $ to transform the 
correlation function of the field $ back to a correlation function of the fields h and e. 
Therefore, we perform the integration 



/ exp 



y^y' 



nd$.d$: 
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Figure 1. CoefScients /^r and of the RGP theory plotted against the chemical 
potential /i for the tunneling rate aJ — k^T. Both functions are symmetric in /i. A 
BEG exists if fijx, > 0. The two points where /ir ~ mark the phase transition between 
the non-condensate and the BEG. The coefficient ^r ist always positive. 



d 



d 



dib;ey) d{hy.el,] 
d d 



exp 



d(hley) d(hy.e\] 



det (3'{1- i') exp 



13" det(l-t') 



[l-i')y,y'+(3'Y.Ke.el,hAl-i') 



y,x' 



exp 



Since we are interested in the limit y — y' ^ oo, and the matrix i^ y, includes nearest- 
neighbour hopping only, the term (1 — does not contribute. This yields for far 
distant lattice points y, y' the expression 

= E {Kexe*M (1 - i'Uyil - i')y,y . 

In this sum, only those terms contribute, where well as x', y' are nearest 

neighbours. In the limit y — y' ^ oo we can assume {blexe*,bx') = (byCyCyfby'Y Thus 
we can use 



E(i-t'k. = E(i-^"'V,x' = (i + i/ 



a) 



We get 



hm (%<!>;,) + hm (bleyel,h 



y-y'->oo 



and therefore 



1 



y-y'^oo 



'y'^y^y'"y 



(1 + y-y'^<x \ y 



(32) 



(33) 



(34) 



On the mean-field level, this justifies our identification of the condensate denity in Eq. 
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Figure 2. Condensate density uq and total density ntot of a translational invariant 
system, plotted against chemical potential for zero-temperature (thick lines) and finite 
temperature (thin lines). 




Figure 3. Condensate density uq = |$p/a^(l + 1/a)^ of a vortex-free condensate in 
a spherical trap with the numerical parameters given in Eq. and different values 
of the chemical potential /i' calculated from the full slave-boson mean-field Eq. ([M]) 
(thick lines) and within RGP approximation (j25p (thin lines). 
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Figure 4. Critical angular velocity plotted against the total number of particles in 
the condensate. 
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Figure 5. Condensate density of a condensate with a single vortex, with same 
parameters as in Fig. [3l calculated from the full mean-field Eq. ([26|l (thick lines) 
and within RGP approximation (|27p (thin lines). The rotating frequencies of the trap 
were chosen to be close to the critical angular frequency Q.c- 




